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Abstract. We present spectroscopy of mid-infrared emission lines in twelve starburst regions, located in eleven 
starburst galaxies, for which a significant number of lines between 2.38 and 45pm were observed with the ISO 
Short Wavelength Spectrometer, with the intention of providing a reference resource for mid-infrared spectra of 
starburst galaxies. The observation apertures were centred on actively star forming regions, including those which 
are inaccessible at optical wavelengths due to high levels of obscuration. We use this data set, which includes 
fine structure and hydrogen recombination lines, to investigate excitation and to derive gas phase abundances of 
neon, argon, and sulphur of the starburst galaxies. The derived Ne abundances span approximately an order of 
magnitude, up to values of ~ 3 times solar. The excitation ratios measured from the Ne and Ar lines correlate well 
with each other (positively) and with abundances (negatively). Both in excitation and abundance, a separation 
of objects with visible Wolf-Rayet features (high excitation, low abundance) is noted from those without (low 
excitation, high abundance). For a given abundance, the starbursts are of relatively lower excitation than a 
comparative sample of HII regions, possibly due to ageing stellar populations. By considering the abundance 
ratios of S with Ne and Ar we find that, in our higher metallicity systems, S is relatively underabundant by a 
factor of ~ 3. We discuss the origin of this deficit and favour depletion of S onto dust grains as a likely explanation. 
This weakness of the mid-infrared fine structure lines of sulphur has ramifications for future infrared missions 
such as SIRTF and Herschel since it indicates that the S lines are less favourable tracers of star formation than 
is suggested by neb ular models which do not consider this effect. 

In a related paper iSturm et al.ll2002l) . we combine our results with spectra of Seyfert galaxies in order to derive 
diagnostic diagrams which can effectively discriminate between the two types of activity in obscured regions on 
the basis of excitation derived from detected mid-infrared lines. 
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1. Introduction 

By definition, starburst galaxies are hosts to sites of re- 
cent star formation with associated rates that cannot be 
maintained over a Hubble time, which can be much higher 
than that determined for the Milky Way. Starburst galax- 
ies are believed to contribute significantl y to the pop- 
ulatio n of massive stars in the Universe ijGallego et alJ 
119951) . Within lOMpc, approximately 25% of high mass 
star formation is attrib uted to only four starburst galaxies 
I Hec kman et al ] |1998Il The reported increase in the star 
formation density o f the Universe for 0.1 ^ z ^ 2 (e.g. 
Madau et ail Il996t iMadau Pozzetti. fc Dickinson! 11991* 
Barger. Cowie. fc Richardsll2000|) and the high frequency 
of starbursts with disturbed morphologies or in interact- 
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ing/merging systems both imply that starbursts are likely 
to play an important role in galaxy formation and evo- 
lution scenarios. By virtue of their proximity, local star- 
bursts can be used as astrophysical laboratories to investi- 
gate the processes that may be ongoing at higher redshifts. 
However, despite intensive observational and theoretical 
modelling campaigns, the properties of local starbursts are 
still not fully understood. In particular, their stellar pop- 
ulations and the state of the interstellar medium remain 
difficult to constrain. 

Observable expressions of a starburst depend upon 
a number of factors, including the stellar initial mass 
function (IMF), star formation history and the age- 
ing/evolution of the stellar population. Both the evolu- 
tion of the stellar population and the spectra of individ- 
ual stars are a function of metallicity, which provides the 
motivation for studying the effects of metallicity that we 
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consider in this paper. A galaxy's interstellar medium af- 
fects nebular diagnostics and, through obs curation, also 
direct stellar ones. The widely used local IScalol l)l986|) 
IMF may not be reflective of the stellar mass content 
for starburst galaxies. IMFs modified at either the low 
or high mass end have been invoked to explain obser- 
vationally derived properties of starbursts such as dy- 
namical masses, spectral line strengths, continuum flux 
and radial velocities (e.g. |Rjeke eLalJll£s5 |Pjjxlevel_alJ 
TflsfllPovon. Joseph, fc Wright 994HDohertv et alllflflfT 



Achtermann k Lacvl Il995t iBeck. Kellv. k Lacvl Il99 



Sternberd Il998t ICoziol. Dovon. k Demers! l200l|) . The 
limit on the upper end of the IMF remains a controver- 
sial issue; mass limits as low as 30 Mq have been proposed. 
However, the presence of high mass stars (up to IOOMq) is 
well known in sites of extreme star formatio n located both 
in th e Galaxy and the local un i verse ( e.g. iKrabbe et al 



19951 ISerabvn. Shupe. k Figel Il998l lEisenhauer et aL 
19981: iMassev k Hunted Il998t iTremonti et alJ l200l|) ~ 



addition, the detection of Wolf-Rayet features in several 
starbursts implies that their prog enitor stars wit h masses 
2; 25 — 60Mq were once present l)Schaererll200ll and ref- 
erences therein). 

Analysis of mid-infrared nebular emission lines shows 
that the excitation of dusty s tarbursts is often low er than 
that of Galactic H II regions l|Thornlev et all EboOl . From 
theoretical mod elling of the [NeIH115 . 5/im / [Nell] 12.8/um 
excitation ratio. iThornlev et alJ (|2000j) suggested that the 
stellar SEDs of 27 starburst galaxies were consistent with 
the formation of massive stars (50 — IOOMq) and argued 
against IMFs with upper mass cut-offs lower than 30A/q. 
The low excitation of galaxies that were originally forming 
high mass stars is then mainly attributed to ageing of the 
stellar population, as opposed to a low upper mass limit to 
the IMF. In addition, the dependence of both stellar evo- 
lutionary tracks and neb ular properties (including excita- 
tion) on metallicity (e.g. lBreso lin 1 Kennicutt, k, Garnettl 



Il999t IThornlev et alJEood: iGiveon et aljl2002a|) suggests 
that elemental abundances must be considered in the anal- 
ysis of the issues detailed above. 

To date, abundance and excitation studies have been 
based primari l y on optical and near-infrared data (e.g. 
01ofssonlll995tlKobulnickv k Skillmanlll996t ICoziol et alJ 



1999; IConsidere et alJl200(]|) . mainly due to observational 
limitations from the ground. Optical abundance studies 
however do not reach the obscured regions dominating 
the activity of many infrared-selected starbursts. Even 
for starbursts which have been well-studied in the opti- 
cal, a combination of radio and infrared measurements 
has conclusively demonstrated that some of the most ac- 
tive star forming sites are optically obscured, with most 
of the bolometric luminosity emerging in the IR (e.g. 
Goriian. Turner, fc Beckl l200ll IVacca. Johnson, k Contil 



2002|) . Thus, investigations of the densest regions of star 
formation are often restricted to the infrared due to ex- 
tinction by dust at shorter wavelengths. 

In general, optical- versus infrared-selected samples of 
starbursts are sensitive to different starburst properties. 



Star formation traced by optical observations is generally 
located within a disk having low to moderate extinction 
with relatively low star formation rates (< 20Mq yr^ 1 ), 
while in the infrared, the densest regions of star forma- 
tion can be probed, occurring as compact events with star 
formation ra tes which may reach <J few x 10 2 Mq yr~ l 
( Riekej200l|). Infrared spectroscopy, probing the dominant 
obscured regions of objects at the dusty end of this se- 
quence, is needed for a complete understanding of the lo- 
cal starburst population with implications for the exten- 
sive populations of starbursts at higher redshift detected 
both in the UV/optical and infrared. 

Here, we present Infrared Space Observatory (ISO) 
Short Wavelength Spectrometer (ISO-SWS) spectra (be- 
tween 2.38 and 45/xm) for 12 starburst regions observed 
in significa ntly more deta i l than most of the galaxies 
studied bv IThornlev et all l)2000|) . a sample with which 
it partially overlaps. This paper complements the nu- 
merous existing abundance analyses in the optical/near- 
infrarcd with the first comparable mid-infrared systematic 
study of a comprehensive sample of starburst galaxies. 
Our analysis is focused upon the gas phase abundances 
of neon, argon and sulphur in this group of local star- 
burst galaxies. These elements emit the strongest emission 
lines in the mid-infrared spectra of starbursts. Deriving 
abundances from a combination of infrared fine structure 
and hydrogen recombination lines has a number of tech- 
nical advantages that minimise uncertainties. Extinction 
in this wavelength range is low, and the derived abun- 
dances (relative to H) are moderately insensitive to varia- 
tions in el ectron temperature and density [abund ances are 
cx T~ - 7 l|Giveon. Morisset. k Sternberd l2002b|) and the 
dependence on n e cancels for abun dances relative to H, see 
Sect. 4.4.1. See also discussion in lForster-Schreiber et alJ 
(200^).]. For all three elements, the major ionisation stages 
found in HII regions are covered by the ISO-S WS spectra, 
which minimises the required ionisation correction factors. 
We relate the results to nebular excitation and thereby 
probe the hot stellar population. 

The data presented include fainter transitions not in- 
cluded in the abundance analysis, and will be useful as 
a reference for work with upcoming mid-infrared spectro- 
graphs on 8m class telescopes as well as observations of 
fainter and/or higher redshift sources with forthcoming 
infrared telescopes such as SIRTF, SOFIA, and Herschel. 

The layout of this paper is as follows: In Sect. 2 we 
describe the observations, the sample selection and data 
reduction. In Sect. 3 we present the line lists and in Sect. 
4 the excitation and calculated abundance analysis. We 
discuss the results in Sect. 5 and finally, our conclusions 
are presented in Sect. 6. 

2. Observations 

2.1. Sample Selection 

Our sample consists of twelve regions in eleven galaxies 
which exhibit starburst characteristics in the infrared. Of 
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the starburst galaxy observations by ISO-SWS present 
in the ISO Data Archive, we selected those for which, 
in addition to fine structure lines, at least one hydrogen 
recombination line was detected. This is a pre-requisite 
since a H-recombination line is used as the reference 
line for the abundance estimation. As mid-infrared H- 
recombination lines are faint in extragalactic sources, this 
requirement therefore also restricts our sample to local 
systems (v m ed = 787fcms _1 , v m ax = 3124/cms -1 ). We 
note that our sample was neither homogeneously selected 
nor is it complete. Therefore, we investigate each galaxy 
on an individual basis and analyse trends found for this 
ensemble of starbursts. Details of the starburst sample are 
given in Table 

2.2. Observations and Data Reduction 

The ISO-SWS instrument on-board ISO performed spec- 
troscopic measurements between 2.38 to 45.2/^m with 
grating spectrometers producing spe ctra of medium 
to hig h resolution \ R=1 00-2000, see Ide Graauw et all 
l|l996|k iLeech et all J2002ft for more details]. Within its 
operational wavelength range lie a number of previously 
seldom-observed emission lines including fine structure 
line and H-recombination lines from which elemental 
abundances may be determined. 

The ISO Data Archive was used to collate all of the 
SWS grating spectra available for each galaxy. Most data 
sets were taken using the SWS02 observation template, 
which provides medium resolution measurements of tar- 
geted lines. In some cases, full grating scans producing 
lower resolution spectra were taken (the SWS01 obser- 
vations template) and a handful of medium resolution 
SWS06 1 measurements were also made. Where multiple 
observations were available, care was taken not to com- 
bine data taken with different aperture centres to ensure 
the same part of the galaxy was being measured for each 
observation. 

We used the SWS interactive analysis (I A3 2 ) to reduce 
the data sets from the edited raw data products. The moti- 
vation for using IA3 to re-reduce the data (rather than the 
pipeline reduction) was to enhance the quality of the re- 
duced product by using interactive tools for dark current 
subtraction, correction between up and down scans and 
outlier masking. The resultant reduce d spectra were then 
subject to further processing in ISAP 3 ijSturm et alll998|) 
to remove further outliers and to perform flat-fielding and 

1 Refer to lLeech et alJ ((2002) for details regarding the ISO- 
SWS and its observation modes. 

2 IA3 is a joint development of the SWS consortium. 
Contributing institutes are SRON, MPE, KUL and the ESA 
Astrophysics Division. A standalone version of IA3 (OSIA) is 
available from http://sws.ster.kuleuven.ac.be/osia/ 

3 The ISO Spectral Analysis Package (ISAP) is a joint 
development by the LWS and SWS Instrument Teams 
and Data Centers. Contributing institutes are CESR, IAS, 
IPA C, MPE, RAL and SRON. ISAP is av ailable from 

:/ /www. ipac.caltech.edu/iso/isap/isap. html 



fringe removal (where necessary) to produce the final spec- 
tral line profiles from which line fluxes were measured. The 
reduced data will also be made publicly available through 
the ISO Data Archive. 

Since the ISO-SWS instrument takes measurements 
in a number of aperture size and detector combinations, 
aperture corrections should be considered for sources 
which are extended with respect to the apertures (of sizes 
14" x 20" , 14" x 27" and 20" x 33" centred on the coordi- 
nates given in Table QJ. As we are interested in emission 
from only the compact nuclear optically obscured infrared 
starburst we used radio, millimetre or infrared measure- 
ments from the literature to determine the starburst ex- 
tent to compare to the ISO-SWS aperture (see Fig. ^ 
for details). We found no sources required correction for 
extension ex cept M82. The corrections for this source are 
described in lFdrster-Schreiber et alJ l|200l|) . 

Error propagat ion through the reduc tion of ISO-SWS 
data i s described in lShipman et al.l (j200 lT> and lLeech et alJ 
(2002). To this error we add an additional 20% calibration 
error to account for systematic inaccuracies in the photo- 
metric calibration (Table 4.7 in lLeech et al-lEoO^ . 

3. Results - The MIR spectroscopic survey of 
local starbursts 

We present a survey of mid-infrared emission lines from a 
sample of infrared starbursts. The emission line fluxes for 
the most commonly detected lines within the range of the 
ISO-SWS instrument may be found in Table F2] 

Our starburst sample encompasses a variety of 
morphological types (Table and includes six 

starbursts exhibiting Wolf - Rayet spectral features 
llSchaerer. Contini. fc Pindaol ^399). From this list 
of Wolf-Rayet ga l axies we exclude NGC 5236 since 
iRosa fc D'Odorica l)l986|) detected Wolf-Rayet signatures 
only in regions centred upon three supernovae sites, none 
of which lie within o u r obse rvation aperture. Recently, 
iBresolin fc Kennicuttl l|2002h report on the possible 
detection of Wolf-Rayet-like features in the spectrum of 
a nuclear 'hot-spot' in NGC 5236. However, due to the 
weakness of the features, for the subsequent analysis, we 
assume the central region does not contain large numbers 
of Wolf-Rayet stars (see also Sect. 15.50 . 

Of our sources, onl y NGC 4945 displa ys evidence for 
harbouring an AGN. ISooon et al.l l)2000|) have demon- 
strated that this AGN does not contribute to the flux 
of the mid-infrared emission lines. In addition, an AGN 
has been advocated to be responsible for the ionisation 
required t o produce radio-recombination lines dete cted in 
NGC 253 l|Mohan. Anantharamaiah. fc Gossll2002jh 

4 This publication makes use of data products from the 
Two Micron All Sky Survey, which is a joint project of the 
University of Massachusetts and the Infrared Processing and 
Analysis Center/California Institute of Technology, funded by 
the National Aeronautics and Space Administration and the 
National Science Foundation 
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Fig. 1. 2MASS quick-look postage stamps (log-scale) of the starbursts with SWS apertures (rectangular) overlaid. 
Additionally, overlaid ellipses indicate the extent of the starburst regions as obtained fro m the literature. Notes on 
indiv idual sources: NGC 253 The MIR emission is dominated by a region of ~10" in size l|Boker. Krabbe. fc Storevl 
1998) and a single super-st ar cluster lying within it jKeto et alll999l) . IC 342 The central radio source is of comparable 
size to the S WS aperture jHummel fc Graevell 19901). II Zw 40 The 12 u m emission originates from a region only 0.5" 
in diameter l(Beck et al.ll2002l) . M82 See lForster-Schreiber et alJ J200lh. NGC 3256 Has a compact double nuclei with 
combined angular size smaller than the SWS aperture l|Norris"feForbedll995|) and is in agreeme nt with the size (<10" 
diameter) determined from MIR mapping at [NeII]12.8//m and the N-band looker et alJll997|) . NGC 3690A & B/C 
The discrete compact r a dio a d CO emitting sources are also the sites of strongest MIR emission - refer to images 
presented in lKeto et al.1 l(l997ll and lSoifer et al.1 (|200lh . The overlaid ellipses are the MIR source sizes d etermined by 
the lat ter authors. NGC 4038/9 The SWS aperture contains emission from SGMC2,3,4,5 as identified bv lWilson et alJ 
( 2000) . NG C 4945 N-band imaging revea ls this source has a compact unreso lved nucleus with surrounding core emission 
oriented along the galaxy's major axis llKrabbe. B5ker. fc Maiolind 12001^ . The overlaid ellipse marks the extent of 
the gas emitting CO(2-l) lines from lOtt et aTT((200lT) . Data from the SWS aperture offset from the radio source 
are not included with in the analysi s in th is paper. NGC 5236 The ellipse marks the starburst extent determined 
from radio ma ps from iTurner fc Hoi \ 19941) (ad ditional hotspots are also mark ed). NGC 5253 extent of [SIV]10.5//m 

r 1^ ^ /l 1 l ^^^^^^^^^^^^^^^^ ^^^^^^^^^^^ l^^^^^s ,\ : /Y rvrrn ml Til attt> 
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HOU, WH 
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act 
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i n 

40 
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1 CO w 

1.02 X 


1U 


TJTT 
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Irr 


NGC 3256 


10 27 51.20 


—43 


54 


13.5 


/\ AAA J 0^7 

0.009487 


36.6 


2.59 x 


1 nil 
10 


HII 


colliding pair 


NGC 3690A 


11 28 33.60 


+58 


33 


46.6 


0.010411 


40. 


3.40 x 
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HII 


pec. merger, WR 


NGC 3690B/C 


11 28 30.90 


+58 


33 


45.6 


0.010447 


40. 


3.42 x 


10 11 


HII 


pec. merger, WR 


NGC 4038/9 E 


12 01 52.90 


-18 


53 


04.0 


0.005477 


21.0 
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HII 


interacting 


NGC 4945 


13 05 27.40 
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05.5 
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4 
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NGC 5236 


13 37 00.50 


-29 


51 


55.3 


0.001721 


5.0 


1.29 x 


10 10 


HII 


SAB(s), face-on, WR 9 


NGC 5253 


13 39 55.80 


-31 


38 


31.0 


0.001348 


4.1 


1.64 x 


10 9 


HII 


Im pec, BCD, WR 


NGC 7552 


23 16 10.80 


-42 


35 


04.2 


0.005287 


21.2 


6.73 x 


10 10 


HII, LINER 


SBbc 



a of SWS pointing centre 

b Luminosity distance calculated assuming Ho = 75fcms~ 1 Mpc~ 1 , Op = 1 

c Infrared luminosity calculated using the presc ription of |Hjatou J _Soife£^ lll985i) 

d WR: Wolf-Rayet features. Further details from lSchaerer. CoTtlni^^Pin(iaoriu^99r fc references therein; BCD: Blue Compact 

Dwarf galaxy; 

e Overlap region containing giant molecular clouds 3, 4 and 5 fro mlWilao n et al.l l|2000F) 
* AGN contribution to mid-IR emission lines negligible llSooon et~gdTl200nli "^ 
9 but refer to Sects. I3lfcl5~5l 



We detect a number of fine structure lines within 
the ISO-SWS spectra. In particular, the fine structure 
transitions of singly and doubly ionised Ar and Ne are 
detected in the majority of sources as well as those of 
[S III] and [S IV] . Other commonly detected transitions are 
[OIV]25.89>m, [FeII]26.0//m and [SiII]34.81/im. The ma- 
jority of these lines are likely to have an origin in elec- 
tron impact excitations in HII regions that have been 
photo-ionised by hot stars (assuming these sources are 
pure starbursts) with a possible contribution from fast ion- 
ising shocks. The [Sill] line is belie ved to partly originate 
in ph otodissociation regions (e.g. ISternberg fc Dalgarnol 
Il995t) at the interfaces between HII regions and the sur- 
rounding dense molecular clouds. The ISO-SWS detected 
lines span a range of excitation potentials from 8.2eV for 
[Sill] to 55eV for [OIV]. The origin of the [OIV] tran- 
sit ion in starburst and Wolf-Rayet galaxies is discu ssed 
in lLutz et all l(l998h and lSchaerer fc Stasiriskal l|l999() re- 
spectively. In addition, the H-recombination lines Brackett 
a, (3 and Pfund a are detected in most of our galaxies. 

We have also measured several transitions of molecu- 
lar hydrogen. H2 data for all the objects in our sampl e 
are included in the analysis of iRigopoulou et all l(2002t . 
Our reductions use the newest calibration files of the 
ISO-SWS instrument and we find the majority of our 
fluxes are consiste nt within the calibration uncertainties of 
Rig opoulou et all and therefore we do not present the Hi 
line fluxes here. We also do not address the broad emis- 
sion and absorption features due to dust and ice that are 
measured in M82 and NGC 253. These are discussed by 
ISturm et all (|2000j) . 



Table 2. Emission line fluxes (uncorrected for extinction), in units of 10 20 Wcm 2 
extinction law and logw < Q(Hq) > calculated as described in the text. 



Also given are the calculated Ay for a mixed model with a 
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NGC4038 


NGC4945 


NGC5236 
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NGC7552 
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3.4 


1.2 


41.0 


4.2 


2.1 


2.1 


0.8 


— 
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2.7 


2.5 


HI Pfd 


3.2970 


2.6 


— 


<4.8 


5.9 


— 


— 


— 


— 


— 


— 


— 


— 


HI Pfg 


3.7410 


3.0 


— 


— 


10.7 


— 


— 


— 


— 


— 


— 


— 


— 


HI Bra 


4.0520 


31.2 


7.0 


2.6 


81.5 


4.4 


2.6 


3.7 


2.5 


9.9 


7.8 


4.5 


4.4 


HI Pfa 


7.4600 


7.6 


2.4 


<4.6 


25.9 


2.3 


1.0 


0.7 


— 


2.3 


— 


2.5 


— 


[MgVII] 


5.5032 


<7.1 


<0.8 


<16.8 


<69.3 


— 


— 


— 


— 


— 


— 


— 


<1.1 


[MgV] 


5.6098 


<3.1 


— 


<17.4 


<79.2 


<2.3 


— 


— 


— 


— 


— 


<1.7 


— 


[Aril] 


6.9852 


184.7 


58.6 


0.2 


267.0 


22.7 


— 


— 


2.2 




— 
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18.5 


[NeVI] 


7.6524 


<0.6 


— 


<6.2 
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<1.3 


<0.4 


<0.5 


— 


<0.7 


— 


— 


— 


[Arlll] 


8.9913 


7.0 


2.6 


2.0 


47.6 


4.7 


1.3 


2.8 


1.8 


<2.3 


— 


5.9 


2.5 


[SIV] 


10.510 


<3.3 


<1.2 


19.8 


14.9 


0.9 


1.0 


3.5 


2.7 


<0.9 


0.9 


34.0 


0.3 


[Nell] 


12.813 


453.6 


90.7 


1.7 


714.0 


89.2 


31.9 


27.7 


7.7 


92.5 


133.9 


7.9 


68.0 


[NeV] 


14.321 


<7.6 


<0.7 


<0.8 


<29.7 


<4.8 


<0.5 


<1.1 


— 


<0.5 


<0.5 


<0.3 


— 


[CM] 


14.367 


7.7 


0.7 


<4.1 


<33.0 


<1.8 


— 


— 


— 


0.3 


1.8 


<1.5 


— 


[Nelll] 


15.555 


32.1 


8.4 


17.0 


126.0 


14.2 


9.3 


20.0 


6.5 


9.8 


6.8 


28.9 


5.1 


[SHI] 


18.713 


99.2 


45.3 


5.1 


252.0 


32.5 


8.2 


18.4 


7.3 


6.3 


54.4 


15.9 


24.6 


[Fell I] 


22.925 


17.1 






46.6 


3.3 








0.7 








[OIV] 


25.890 


6.7 


0.7 


1.1 


8.0 


0.8 


0.6 


1.0 


<0.3 


3.5 


0.8 


0.9 


1.0 


[Fell] 


25.988 


14.2 


4.8 


1.0 


36.7 


2.8 


1.6 


2.2 


0.3 


3.5 


4.2 


0.9 


2.0 


[SIII] 


33.481 


197.3 


89.3 


11.4 


562.0 


56.4 


27.1 


34.1 


16.0 


51.4 


109.8 


34.4 


41.1 


[Sill] 


34.815 


273.2 


131.6 


4.5 


839.0 


84.0 


43.4 


42.2 


10.5 


107.4 


138.0 


16.7 


61.5 


[NcIII] 


36.013 


<4.5 


<2.3 


1.4 


18.8 


2.4 


<1.3 


1.8 


1.8 


<1.6 


1.0 


4.9 


<0.8 






8.67 


7.22 


10.36 


52.00 E 


<8.73 


38.94 


<57.13 


37.34 


188.19 


5.00> 


5.08 


<73.90 


log w < Q{H ) > 




53.02 


51.94 


53.77 


53.83 


54.53 


54.62 


54.76 


53.81 


53.61 


52.90 


52.34 


54.40 



a Additional Lines 
6 Additional Lines 



xlO 
xlO 

Additional Lines (xlO 



- 20 Wcm 
- 20 Wcm 
- 20 Wcm 



[FeII]5.34017Atm 32.4, [PIII]17.8850/im 3.0, [FeII]17.9359/mi 9.4 

[FeII]5.34017/im 43.2, [PIII]17.8850/im 18.0, [FeII]17.9359/xm 5.9, [FeIII]33.0384/im < 20. 
[FeII]5.34017/im 5.0, [PIII]17.8850/im < 0.5, [Fell] 17.9359/^1 0.7 



d This [Aril] flux is measured from an SWS AOT01 spectrum which is not centred on the radio starburst, therefore Ar fluxes and abundances for NGC 4945 are not used 
the subs equent analysis. 



from 
from 



Forster-Schreiber et all <|200ll) 
Genzel et alJ ljl998h 
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4. Analysis 

4.1. Extinction Correction 

Several of our starburst galaxies are known to be dusty 
systems with large obscuration; therefore the observed 
line fluxes must be corrected for extinction. To deter- 
mine the correction, we used the ratios of Br (3 /Bra 
and [SI II] 18.7/33.5. The intrinsic emissivities of H- 
recombination lines are theoretically well determined and 
may be compared to the ISO-SWS measurements to de- 
termine the level of extinction. Moreover, ratios of H- 
recombination line emissivities are relatively insensitive 
to variations in electron temperature and density and 
therefore reasonable (rather than precisely accurate) es- 
timates of the two quantities is sufficient to identify 
the a ppropriate intrinsic line ra tio. Emissivity coefficients 
from lStorev fc Hummed 1|1995|) for case B recombination 
with an electron temperature of 10 4 K and density of 
100cm -3 were used to calculate the intrinsic line ratios 
for the H recombination lines. We chose no t to use the 
H-recombination line ratio of Bra/ P fa since lLut3 ( 19991) 
suggest that the extinction at 7.5/rni may be higher tha n 
predicted by standard extinction curves (|DraineJ [1989), 
with the possibility of significant variation among sources. 

The ratio [S'///]18.7/33.5 is more sensitive to density 
than the H recombination lines but there are good reasons 
to assume that the starbursts are close to the low density 
limit (Sect. 14.2. 1[) . We solved the rate e quations for this 
ratio using the updated atomic data of iTaval fc Guptal 
( 199^) and confirmed that the resulting ratio is relatively 
insensitive to temperature. We adopt a value of 0.5 for 
densities close to the low density limit. 

The extinction was calculated based upon a model 
in which obscuring material and emitting gas are uni- 
formly mixed. Such a 'mixed' model is considered more 
appropriate for a starburst region observed with our large 
apert ures than a simple foreground screen ijMcLeod et alJ 
Il993|) . The attenuation of the emitted intensity may be 
expressed as 



1 



(1) 



where t\ is the optical depth at the wavelength of in- 
terest (A) which is proportional to the optical depth in 
the V band Ty through the extincti on law. We assumed 
the Galactic centre exti nction law of | Lut d ljl999l) between 
3-8.5^im and adopted a iDrainei (ll989F ~ extinction law at 
all other wavelengths, with a normalisation of the silicate 
feature of tq.q/tv — 0.122. When comparing two recombi- 
nation lines, the observed ratio may be expressed in terms 
of ry using the appropriate emissivity and extinction val- 
ues. 

In the case that all lines of both ratios Brj3/Bra 
and [S7i7]18.7/33.5 were detected, we quote a weighted 
average of the calculated Ay. For sources where only 
one of these ratios was measured, the corresponding Ay 
was adopted. For NGC3256, NGC3690B/C and NGC7552 



both measured ratios were higher than the intrinsic ratios, 
implying Ay — which is unrealistic for our dusty star- 
bursts. We attribute the greater than intrinsic ratios to the 
large errors on our flux measurements (~ 20%). Therefore, 
as no meaningful values for the extinction could be found 
for these sources, we replaced the measured ratio by the 
measured ratio minus its associated error to determine a 
one sigma limit on Ay. For these sources we apply an 
extinction correction using the one sigma limit of Ay on 
the data but indicate the location of the uncorrected data 
(Ay = 0) in the figures presented in the following sections. 
The estimated Ay used to correct all line fluxes are given 
in Table H 

4.2. Nebular Conditions 

4.2.1. Electron density 

To determine the electron density of our sources 
we used the ratios of [SIII]18.7/^m/33.5/im and 
[NeIII]15.6/^m/36/im and compared them to theoret- 
ical values obtained from a solution to the rate equations, 
obtaining relations of line ratios a n d de nsities close 
to those shown by I Alexander et alJ l)l999l) . partly on 
the basis of earlier atomic data. Assuming a one-zone 
ionisation model, we found that all of our starbursts 
lie in the low density limit (i.e. well below the critical 
density) within the errors. Using only [SIII] both in 
the context of density indication and in extinction 
correction estimates (Sect. 14 - If) allows degenerate low 
extinction/low density and high extinction/high density 
solutions. Low densities found from [Nelll] and gen- 
eral consistency of the recombination-line based and 
[SIII]-based extinction estimates assure however, that 
it is reasonable to adopt low density conditions for the 
entire sample. Our computations of collisional excitation 
indicate n e ~ 10 — 600cm -3 therefore we assumed an 
average electron density of 300cm -3 . We were then able 
to use the expressions in the following section to calculate 
the nebular abundances in the low density limit (within 
which collisional de-excitation may be ignored). 

4.2.2. Electron Temperature 

As with density, electron temperature affects emissivi- 
ties and thus inferred abundances albeit with only mod- 
est dependences for fine structure transitions and H- 
recombination lines. Since abundances are calculated as 
ratios of fine structure line to H-recombination line emis- 
sivities the weak depend ence on electron temperature 
partly cancels out (e.g., iGiveon. Morisset. fc Sternberg! 
2002bj). Accurate determinations of electron tempera- 
ture can be obtained from recombination lines at ra- 
dio and millimetre wavelengths. Measurements in the 
latter wavelength range are easier to interpret than in 
the radio since stimulated emission effects (which are 
present in the radio) are negligible and the millime- 
tre line flux primarily originates in spontaneous emis- 
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sion. Therefore we sought electron temperature determi- 
nations based upon millimetre observations from the lit- 
erature; for the starburst galaxies M82 and NGC 253, 
electron temperatu res of T P ~ 5000 if zfa 1000 j^T have 
been d etermined bv lPuxlev et alJ l|l989f) and lPuxlev et al.1 
l)l997(l . An electron temperature estimate based upon ra- 
dio recombination lines is also availab le for NGC 5253 
ijMohan. Anantharamaiah. fc Gossll200l[) . as well as esti- 
mates derived fro m optical spectroscopy for NGC 5253 
and II Zw 40 llCamobell. Terleyich. fc Melmckl [li 



IVacca fc Cont3ll992t IWalsh fc Rovlll993h The results for 
these two objects agree well with the radio giving electron 
temperatures - 10000 - 12500if. 

The known electron temperatures in our starburst re- 
gions span a factor ^2 consistent with the well known 
corre lation between electron temperature a nd metallicity 
(e.g. ICampbell. Terlevich. fc Melnickl Il98l for a sample 
of on-average lower metallicity objects). For the remain- 
ing starburst galaxies, no accurate radio- or millimetre- 
based temperature determinations are available. We have 
chosen not to rely on electron temperature estimates de- 
rived from optical spectroscopy since these may not reflect 
electron temperatures of our obscured infrared-radio star- 
bursts which cannot be probed by optical spectroscopy. 
For this reason, and to perform inter-comparisons between 
the sources without introducing a bias due to uncertain 
adopted electron temperatures, we elected to use a sin- 
gle 'representative' temperature for the entire sample of 
galaxies. Since M82 and NGC 253 are archetypal star- 
bursts we chose our 'representative' electron temperature 
to be Tf- = 5000AT. However, in doing so we appreciate 
that deviations from this 'representative' temperature are 
present and will affect absolute abundances, thus distort- 
ing the trends derived in the inter-comparisons. We will 
indicate these effects (which are largest for high T e ob- 
jects) in the appropriate location. In particular for the 
BCD galaxies, the magnitude of the effect of using T e of 
~ lOOOOiiT on the abundances is shown in the relevant 
plots. 

For consistency we have also adopted the same 'repre- 
sentative' el ectron temperature (T j 1 = 5000-fQ for the H II 
regions from lGiveon et all l|2002al) , which we use as a com- 
parison sample to the starbursts. We have re-calculated 
extinctions and abundances for this data set in order to 
eliminate any discrepancies due to different assumptions. 
The H II regions and starbursts are treated differently only 
in the obscuring model used to determine the extinction 
correction. For HII regions, often located in the Galactic 
plane but far from the Sun, we have adopted a uniform 
foreground screen model since we assume that extinction 
due to the foreground galactic ISM exceeds the extinc- 
tion due to dust mixed within the ionised gas of the H II 
region itself. Similarly to the BCDs, the assumed repre- 
sentative temperature (5000K) is below the values tha t 
have b een previously used [10000K llGiveon et alJl2002a|) . 
7500K l)Martm-Hernandez et al.l2002a|) ]. Therefore the ef- 
fect of increasing the temperature to 10000K shown on 
the figures in the next section is also applicable to part of 



Starbursts 
• WR Starbursts 
HII R Col < 7.2kpc 
HII R Gol > 7.2kpc 
Extragalactic HI 



<■ c l <yf C525j 

C4038 1 



sB 



1 .00 

[Nelll]15.5/[Nell]12.8 



Fig. 2. Excitation diagnostic: [Nell I] 15.5/ [Nell] 12.8 
versus [ArIII]8.9/[ArII]6.9. Starbursts are plotted as 
filled squares with the Wolf-Rayet starbursts marked by 
an enclosing diamond. Galactic HII region data from 
iGiveon et alJ fc002al) are separated at the median galac- 
tocentric radius to provide 'inner' (open triangles) and 
'outer' (open circles) samples. Finally, we plot the extra- 
galactic H II regions from the same reference as open stars. 
The same key is used to distinguish the data points in all 
subsequent plots. Within this diagram low excitation is 
to the bottom-left quadrant and high excitation is to the 
top-right quadrant. All fluxes involved in this diagram 
were corrected for extinction. For NGC3256, NGC3690 
B/C and NGC 7552 the effect of adopting A v = on the 
derived quantities is plotted as an open square. 



the IGiveon et al.1 l|2002a[) sample of HII regions, namely 
those at large Galactocentric radii or inside the Magellanic 
Clouds. 

4.2.3. Lyman continuum photon emission rate 

The Lyman continuum photon emission rate [Q(Hq)] is 
the number of photons emerging per second from a star 
cluster with sufficient energies to ionise hydrogen. We have 
used H recombination lines to determine the average hy- 
drogen ionising rate [Q(Hq)] using 



Q(H ) 



a B (T e ) 



hvu 



(Te) 



(2) 



The ratio of recombination coefficients {&b ) is 
relatively insensit ive to temperature (see equation 5.23 
lOsterbrockl Il989 j ). The recombinat ion coefficients were 



taken from lStorev fc Hiunmerl (1995). The calculated val- 
ues of Q(H ) are given in Table [3 

4.3. Excitation 

We may investigate excitation within the starburst sam- 
ple using ratios of fine structure line fluxes of different 
ionic species of the same element (e.g. X +l+1 / ' X +t ) . This 
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10.00 r 




0.01 0.10 1.00 10.00 100.00 

[Arlll]15.5/[Arll]12.8 

Fig. 3. Excitation diagnostic: [ArIII]8.9/[ArII}6.9 ver- 
sus [SIV] 10.5/ [SIII] 18.7. Key as in FigureEl Within this 
diagram low excitation is to the bottom-left quadrant and 
high excitation is to the top-right quadrant. All fluxes in- 
volved in this diagram were corrected for extinction. 



10.00 r 




0.01 0.10 1.00 10.00 100.00 

[Nelll]15.5/[Nell]12.8 

Fig. 4. Excitation diagnostic: [NeIII]15.5/[NeII]12.8 
versus [SIV] 10.5/ [SIII] 18.7. Key as in Figure^ Within 
this diagram low excitation is to the bottom-left quadrant 
and high excitation is to the top-right quadrant. All fluxes 
involved in this diagram were corrected for extinction. 



We have formed diagnostic excitation planes by com- 
bining the fine structure line ratios as shown in Figs. EJ 
13 and 0] The starbursts are relatively spread out across 
all three excitation planes. The majority of sources oc- 
cupy the lower left corner of the planes, which indicates 
that these have low excitation. Extending to higher exci- 
tation in both axes are the BCD galaxies and those that 
exhibit Wolf-Rayet features (Wolf-Rayet galaxies are dis- 
cussed in more detail in Sect. 15.5]) . This result is consis- 
tent with existing knowledge of bot h of these source types 
l|Schaerer. Contini. fc Pindaolll999|) . 

Also plotted in the excitation plane s are extinc- 
tion co rrected data of HII regions from iGiveon et alJ 
il2002al) . The authors found an increasing excita- 
tion gradient with galactocentric radius within 
the galaxy (see also [Martin- Hernandez et alJ l2002at 
IGiveon. Morisset. fc Sternberg! l2002bh . Therefore these 
data have been separated into two broad bins based upon 
galactocentric radius (denoted by 'inner' and 'outer') and 
reflect low and high excitation. The distribution of the 
starbursts is well matched with the H II region data. The 
locations of the 'inner' HII regions is similar to those 
of low excitation starburst galaxies. The exceptions are 
the Wolf-Rayet and BCD galaxies which are located in 
regions similar to those of the 'outer' and extragalactic 
HII regions. These extragalactic HII regions have been 
heterogeneous l y selec ted from the Magellanic Clouds 
I Giveon et ail 120024 - Th e SMC and LMC are both 
low metallicity dwarf galaxies containing Wolf-Rayet 
stars and thus the coincidence of our Wolf-Rayet and 
BCD starbursts with the extragalactic HII regions is as 
expected. 

For both the starbursts and the HII regions we see a 
general correlation between the three excitation ratios of 
Ne, Ar and S in Figs. 121 II The existence of a correlation 
between all three ratios implies that they probe the same 
excitation mechanism. Moreover, the ratios are sensitive 
to photons emitted 2, ^ 2.5 and ^ 3 Rydbergs and 
thus the correlations suggest that regions with ionisation 
potentials between 2 — 3 Rydbergs in the ionising spectra, 
are likely to be closely related in 'hardness'. 



ratio, to the first order and for a given ionisation parame- 
ter, is proportional to the number of photons which have 
sufficient energy to ionise the state X +t relative to the 
number of Lyman continuum photons. The difference in 
ionisation potentials between the ionic species means that 
ratios of these lines are sensitive to the shape of the spec- 
tral energy distribution (SED) of the ionising source. We 
have used (and combined) ratios of some of the strongest 
lines, [NeIII]15.5/[NeII}12.8, [ArIII]8.9/[ArII}6.9 and 
10.5/ [SIII] 18.7, from our mid-infrared survey. Ar + , 
S ++ and Ne + have ionisation potentials corresponding to 
~ 2, ~ 2.5 and ~ 3 Rydbergs, respectively. Therefore 
these line ratios can constrain the ionising spectra over a 
range of a few Rydbergs. 



4.4. Elemental Abundances 

In the following, we derive the elemental abundances for 
our sample galaxies and discuss their relation to other fac- 
tors, in particular the role of metallicity in determining the 
wide range of excitations. Although we detect fine struc- 
ture lines for several elements in some objects (Table EJ, 
we restrict our abundance analysis to only those lines orig- 
inating from H II regions and detected in the majority of 
sources within our sample i.e. neon, argon and sulphur, 
which are all thought to be 'primary' products of nucle- 
osynthesis. They are products of oxygen burning, possibly 
in th e late stages of evolution of massive stars (|Garnettl 
2002). Thus the abundance of these primary products are 
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likely to be correlated with the abundance of oxygen and 
should trace metallicity well. 

4.4.1. Prescription 

We calculate abundances using the following prescription. 
In general, the flux (F x ) of a given ionic species (X +1 ) 
may be written as the product of the number density of 
the ionic species (n x ) , the electron density of the ionised 
gas (n e ) and j x - the emissivity of the species, i.e., 



F H + 



x + i ,Y+i 



(3) 



The ionic abundance (n x+ /n H+ ) may then be cal- 
culated from the flux ratio of the fine structure lines of 
an ionic species to a H-recombination line, multiplied by 
the ratio of the emissivities. In HII regions nearly all the 
hydrogen is ionised, therefore n H+ ~ n H . The emissivity 
of the H-recombinati o n line is obtained from the tables 
of IStorev fc Hummerl (jl995h . Strictly speaking, the emis- 
sivities of ionic fine structure lines should be obtained by 
solving the rate equations for a full multilevel system un- 
der the appropriate conditions. However, for ground state 
fine structure transitions and a low density H II region, the 
following simple formula provides a good approximation 



E 



he 

T 



8.629 x 10" 



-6\ 1/2 



T, 



1/2 



tot 



kT l: 



(4) 



where j has units of erg cm 3 s _1 , T; u is the effec- 
tive collision strength, loi is the statistical weight of the 
lower level, \iu is the energy difference to the collision- 
ally excited level from the ground state, and the summa- 
tion is over those fine structure levels which, when col- 
lisionally excited from the ground state, decay through 
the line of interest. For the lines we used (Table |3J), it 
is sufficient to consider one upper level only, with the ex- 
ception of the [NeIII]15.5/im and [ArIII]9.0/im lines where 
two levels are needed. Effective collision strengths were ob- 
tained from the appro priate reference of the IRON project 
( Hu mmer et al 1ll993t Table©. The abundance of an ionic 
species may be calculated using 



F 



x+i 



F H + 



(5) 



Once an ionic abundance has been estimated using 
the formula above (see Table 0J, we sum the abun- 
dances for each ion from a given element and apply 
the ionisation correctio n factor s (hereafter ICF) from 
IMartm-Hernandez et alJ l|2002a[) to correct for the con- 
tribution from u nobserved ionic species. The ap plication 
of the ICFs from IMartm-Hernandez et alJ l|2002a[) is ap- 
propriate for these data since we are using the same ionic 
species to determine the abundance, and the model grid 
used to calculate the ICFs is sufficient to include the 



10Q00K 5Q00K 
< 




[Ne/H] 



Fig. 5. Excitation against abundance for neon. Starbursts 
are plotted as asterisks with the Wolf-Rayet starbursts 
marked by an enclosing diamond . Galactic and local- 
extragalactic HII region data from iGiveon et a"fl j.2002a) 
are also plotted as open triangles, circles and stars. For 
a given metallicity the starbursts are of lower excitation 
than the HII regions. The arrow indicates the effect on the 
abundances of changing the adopted electron temperature 
from 5000 to 10000 K. 



sources similar to ours (i.e. those in the low density limit). 
Since the two dominant ionisation stages of Ne, Ar and 
S for HII region conditions are observed in our galax- 
ies, the ionisatio n corrections factors are always close to 
unity. Following iMartm- Hernandez et al.1 I 2002aj) . we do 
not correct our Ne a bundances as the ICF is ~ 1 (see also 
IGiveon et al]l2 002a) and for S we applied a factor of 1.15. 
The ICF for Ar is seen to increase with excitation, from 
1 to a maxi mum of 1.35. We interpolate over the factors 
calculated bv lMartm- Hernandez et all 1 2002ai Fig. 12) to 
obtain an ICF appropriate for the excitation (as deter- 
mined from the Ne ratio) of each of our sources. 



4.4.2. The abundance of neon and argon 

In Figs. El and HJI we show the relationship between 
excitation and abundance 5 for Ne and Ar for the star- 
burst sample 6 . A trend of decreasing excitation with 
increasing abundance is clearly seen in Ne and Ar for 



5 We calculate our abundances relative to solar 
iGrevesse fc Sauvall ll998T> to be consistent with extra- 
galactic abundance work. It is widely thought that the Sun 
is under-abundant i n heavy metals co mpared to its local 
neighbourhood (e.g., ISnow fc WitdlT996t) . Adopting different 
reference abundances will only change the absolute values but 
not affect the correlations discussed in the following. For infor- 
mation, we show in Figs. l71lUl the effects of different reference 
abundances by arrows marke d 'ISM' (inters tellar medium 
[ Wilms. Allen, fc McCravl l2000t) and 'Orion' JSimpson et alJ 

6 For the starburst regions A and B/C in NGC 3690, argon 
abundances are not presented since no measurement exists for 
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Table 3. Atomic data used to calculate the fine structure line emissivities and abundances. Effective Collision Strengths 
are interpolated to a temperature of 5000K. In addition the Bra r ecombination-line was used, adopting an emissivity 
of ef ™ = 2.298 x 1Q- 26 erg cm 3 s" 1 from lStorev fc Hummer! l|l995[) for n e ~ 100cm"3 and T e ~ 5000K. 
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[Nell] 


12.8136 
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4 
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SaraDh fc Tullv (1994) 


[Nelll] 


15.5551 


0.730 


5 


0.07971 


Butler & ZeiDDen ( 1994) 






0.199 
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0.11413 


Butler fc Zeippen (1994) 


[Aril] 


6.9853 


2.70 


4 


0.17749 


Pelan & BcrrinEton (1995) 


[Arlll] 


8.9914 


3.182 


5 


0.13789 


Galavis, Mcndoza, & Zeippen (1995) 






0.663 


5 


0.19469 


Galavis, Mcndoza, & Zeippen (1995) 


[SIII] 


18.7130 


1.41 


1 


0.10329 


Taval & Gupta (1999) 


[SIV] 


10.5105 


8.10 


2 


0.11796 


Saraph & Storcv (1999) 



Table 4. Calculated elemental abundances for each source, per ionic sp ecies and the total (ICF corrected) abun- 
dance with reference to solar abundances ([X/H]): for Ne, Ar and S from iGrevesse fc Sauvall lll998l) (first line). The 
abundances presented are for T e = 5QQQK. See text for the effect of this assumption on the lowest abundances. 
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2.4 ±0.7 


12.1 


± 


3.8 


6.6 


± 


2.6 


5.1 ± 


1.5 


9.4 


± 


3.0 




6.8 ±2.1 


II Zw 40 


1.4 ±0.8 


59.0 ±20.8 


0.6 ±0.2 


0.1 


± 


0.1 


14.6 


± 


6.0 


0.9 ± 


0.3 


2.9 


± 


1.1 


32.3 ± 10.8 


4.4 ± 1.1 


M82 


17.0 ± 7.2 


11.6 ±4.9 


1.5 ±0.6 


4.4 


± 


1.9 


15.8 


± 


6.7 


2.4 ± 


0.8 


4.3 


± 


1.8 


1.2 ± 0.5 


3.2 ± 1.3 


NGC253 


31.0 ± 9.4 


9.3 ±3.1 


2.7 ±0.8 


8.6 


± 


2.7 


4.1 


± 


1.5 


3.6 ± 


1.1 


4.6 


± 


1.5 




3.4 ± 1.1 


NGC3256 


43.2 ± 18.5 


29.1 ± 13.1 


3.8 ± 1.5 


7.4 


± 


3.3 


19.3 


± 


8.9 


3.8 ± 


1.4 


10.7 


± 


4.7 


0.8 ±0.5 


7.8 ±3.4 


NGC3690A 


23.8 ± 10.1 


27.6 ± 11.7 


2.2 ±0.8 








12.9 


± 


5.5 


0.5 ± 


0.2 


4.4 


± 


1.9 


2.4 ± 1.0 


3.4 ± 1.4 


NGC3690B/C 


14.4 ±6.1 


40.2 ± 17.0 


1.5 ±0.5 








21.1 


± 


9.0 


1.0 ± 


0.4 


7.0 


± 


3.0 


6.4 ±2.7 


5.5 ±2.2 


NGC4038 


6.1 ± 2.6 


20.3 ±8.6 


0.7 ±0.2 


1.2 


± 


0.5 


17.6 


± 


7.5 


1.4 ± 


0.4 


4.2 


± 


1.8 


6.6 ±2.8 


3.5 ± 1.3 


NGC4945 


17.2 ±6.5 


6.7 ±2.6 


1.5 ±0.5 


















0.9 


± 


0.4 




0.6 ±0.3 


NGC5236 


36.8 ± 15.6 


8.0 ±3.4 


3.1 ± 1.3 


















10.1 


± 


4.3 


0.4 ±0.2 


7.4 ±3.1 


NGC5253 


3.8 ± 1.7 


59.4 ±22.3 


0.8 ±0.2 


0.7 


± 


0.4 


22.2 


± 


9.0 


1.5 ± 


0.5 


5.2 


± 


2.0 


28.1 ± 10.5 


5.8 ± 1.6 


NGC7552 


29.5 ± 12.5 


8.4 ±3.6 


2.5 ± 1.0 


5.6 


± 


2.4 


16.3 


± 


6.9 


2.9 ± 


1.0 


7.9 


± 


3.3 


0.4 ±0.2 


5.7 ± 2.4 



I10000K 5000K 




[Ar/H] 



Fig. 6. Excitation against abundance for argon. Symbols 
are as in Fig. |3] 



the starburst regions. The Wolf-Rayet galaxies have 
the highest excitation and the lowest abundance. In 



the strong [Ar 1 1]6. 99 /ira line. Therefore any calculated abun- 
dances for these regions would be severe underestimations. 



general, the non- Wolf-Rayet starbursts in our sample 
have super-solar metallicities (X/H) ~ 1 — 3(X/H)q 
(where X is Ne or Ar) whereas the Wolf-Rayet and BCD 
galaxies are sub-solar. As discussed in Sect. 14.2.21 our 
adopted 'representative' electron temperature of 5000 K 
is a good assumption for high metallicity starbursts 
and high metallicity HII regions. However, the same 
is not true for the BCD galaxies and the lower metal- 
licity H II regions among the sample of iGiveon et alJ 
(2002a). Therefore, the effect of changing the electron 
temperature from 5000 to 10000 K on the calculated 
abundances is shown with an arrow in the diagrams. 
By adopting a more appropriate electron temperature, 
the absolute infrared-derived abundances of the low 
metallicity BCDs are close to the optically-derived values 
[0.2(X / ff^l (e.g. ISunter. Gallagher, fc Rautenkranj 



1983 IVacca fc Contil ll 992t IWalsh fc Rovl 

Masegosa. Moles, fc Campos-Aguilarl Il994l) . Overall, 



the effect on the excitation/abundance diagram of Fig. 
[S] of adopting individual electron temperatures (which 
are unavailable for the entire sample) for each starburst 
would be to tilt the correlation: the high excitation / low 
metallicity end would move to ^2 times lower abundance, 
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Table 5. Neon abundances for our sample of galaxies cal- 
culated assuming a fixed electron temperature of 5000K 
(column 2) and individual electron temperatures (3rd col- 
umn). For the latter, abundances are derived from inter- 
polating between objects with known electron tempera- 
tures (values indicated by bold face) , on the basis of their 
[NeIII]15.5/[NeII]12.8 excitation ratios. 



Name [Ne/ H](5000K) [Ne/H](T e ) 



IC342 2.4 2.1 

II Zw 40 0.6 0.3 

M82 1.5 1.5 

NGC253 2.7 2.7 

NGC3256 3.8 2.1 

NGC3690A 2.2 2.0 

NGC3690B/C 1.5 1.9 

NGC4038 0.7 1.9 

NGC4945 1.5 2.1 

NGC5236 3.1 2.1 

NGC5253 0.8 0.4 

NGC7552 2.5 2.1 




1 



[S/H] 

Fig. 8. Excitation against abundance for Sulphur. 
Symbols are as in Fig.0 Unlike the corresponding figures 
for neon and argon, no correlation is seen. 



while the high metallicity/low excitation end remains 
unchanged. 

From the literature we found accurate electron tem- 
peratures for the infrared starburst regions of M82, NGC 
253, II Zw 40 and NGC 5253 (as described in Sect. l4~2~2l 
which we used to re-calculate the elemental abundance 
of neon (Table EJ in these galaxies. II Zw 40 and NGC 
5253 (T e ~ 12500K and - llOOURT, respectively) lie at 
the low abundance, high excitation end of our abundance- 
excitation correlation, while M82 and NGC 253 (T e ~ 
5000-ftT) lie at the opposite extremity. We then defined a 
linear relationship between excitation and true Ne abun- 
dance (calculated using the measured actual electron tem- 
peratures) defined by these two pairs of galaxies. We could 
then estimate neon abundances for the remaining galaxies 
in the sample by interpolating along this relation based 
upon their [NeIII]15.5/[NeII]12.8 excitation ratios. This 
process yielded abundance estimates that are essentially 
independent of our knowledge (or lack thereof) the true 
electron temperature for each source. The estimated abun- 
dances ( Table for the low excitation sources are all close 
to - 2Z Q . 

As 'primary' products of stellar nucleosynthesis, the 
abundances of neon and argon are expected to be corre- 
lated. The overall abundance relation between Ne and Ar 
is shown in Fig. which indeed shows a clear correlation 
between the two elements that holds for all objects in our 
sample, from the BCDs and WR galaxies to the dusty 
starbursts. This implies that the ISM of these galaxies 
has been enriched by similar processes for both elements. 
The [Ar/Ne] ratio shown in the right panel (for which the 
T e effects for the two elements largely cancel) suggests an 
above solar average argon to neon abundance ratio with 
no obvious trend among the galaxies. 



4.4.3. The abundance of sulphur 

We discuss the abundance of sulphur separately since we 
do not see the clear correlation between excitation and 
abundance that is seen in Ne and Ar for the starbursts 
(Fig. [SJ. In addition, all galaxies appear to have sub-solar 
sulphur abundances, with NGC 4945 being anomalously 
low. Furthermore, on comparison to the abundances of Ne 
and Ar, we find that the galaxies not showing Wolf-Rayet 
features have lower S abundances (< [S/H] >= 0.6) than 
expected from their typically high Ne and Ar abundances 
(< [Ne/H] >= 2.4, < [Ar/H] >= 2.7). Yet the S abun- 
dances for the Wolf-Rayet galaxies (which typically have 
lower Ne and Ar abundances) are similar to that expected 
from their Ne and Ar abundances {[X/H] ~ 0.2). 

We confirm this result in Figs.l^land llOl where we in- 
vestigate the correlations between S abundances and Ne 
and Ar. S is also expected to be correlated with Ne and 
Ar since they are all primary products of stellar nucle- 
osynthesis. However, we see no correlation between the 
abundances of S in the galaxies with either Ne or Ar. The 
right panels of Figs. 151 and 1101 compare the abundance ra- 
tios [S/Ne] and [S/Ar] with the Ne and Ar abundance. 
They reflect the trend towards a stronger sulphur 'defi- 
ciency' at increased Ne or Ar abundance noted above. S 
is of ten referred to as an excellent tracer of metallicity 
(e.g. lOev fc Shields! l2000j) and in studies of the ISM is 
commonly used as a non- or slightly-deplet ed (^ 30%) el- 
ement (e.g. lWilms. Allen, fc McCravlEoob]) . However, the 
significant underabundance of S implied by the results of 
this study clearly has implications for the use and inter- 
pretation of sulphur fine structure lines and sulphur abun- 
dances in starburst sources. This result suggests that S is 
not a good tracer of metallicity in starburst galaxies. We 
discuss possible explanations for this depletion of sulphur 
below. 
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[Ne/H] 



[Ne/H] 



Fig. 7. Left: [Ne/H] abundance vs. [Ar/H] abundance. The key is as in Fig. El Right: [Ar/Ne] abundance ratio vs. 
[Ne/H] abundance. The key is as in Fig. El In both panels, a dashed line represents the solar abundances and abundance 
ratios for these elements. The arrows marked 'ISM' and 'Orion ' denote the effect on the data of using different reference 
abun dances other than sola r - 'ISM' interstellar medium from lWilms. Allen, fc McCravl 1 200(t l and Orion abundances 
from ISimoson et all l)l998l) . Adopting these different reference abundances only change the absolute values but not 
the correlation. 




Fig. 9. Left: [Ne/H] vs. [S/H] abundance. The key is as in Fig.El Right: [S/Ne] v [Ne/H] abundance. The key is as in 
Fig. El For both panels, dashed lines represents the solar abundance and abundance ratios. 




[Ar/H] 




[Ar/H] 



Fig. 10. Left: [Ar/H] v [S/H] abundance. The key is as in Fig. El Right: [S/Ar] v [Ar/H] abundance. The key is as in 
Fig. El For both panels, dashed lines represents the solar abundance and abundance ratios. 
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5. Discussion 

5.1. Comparison to optical/NIR derived abundances 

A direct comparison to the abundances derived from opti- 
cal/NIR studies is far from straightforward. In particular, 
our obscured starbursts are not directly comparable to 
optical derived abundances since the regions we are prob- 
ing may be completely obscured. Moreover, when making 
comparisons for our sources of large angular size one must 
ensure that the spectra in both regimes are probing the 
same regions. 

The abundances obtained from our infrared study 
complement optical/NIR determinations by adding Ne, 
Ar and S which have predominant excitation species in 
the mid-infrared wavelength range. We searched the lit- 
erature for abundances determined from optical spec- 
troscopy for these elements, or for oxygen with which 
they are expected to be correlated. For II Zw 40 
and NGC 5253, optical based abundances are con- 
sistent with our mid- infrared results (Z ~ 0.2 Zq 
iKobulnickv et alJ Il997|) . For our low excitation star- 
bursts with metallicities close to 2Zq, while optically de- 
rived values are in agreement that the starbursts have 
super-solar metallicities, the estimates are slightly dis- 
crepant with our infrared de terminations. For NGC 3256 
l|Bresolin fc Kennicuttl 12 002. and references therein) and 
NGC 7552 l|Storchi-Bergman. Calzetti & Kinnevl Il994f) 
the derived metallicities are slightly lower Z ~ 1.5Zq, 
whereas for N5236 the determinations is slightly high er 
Z ~ 2.6Z l|Storchi-Bergman. Calzetti fc Kinnevlll994^ . 

5.2. Excitation and Abundance 

The strong negative correlations seen for Ne and Ar in 
Sect. I4.4?2*l indicate that 'spectral hardness' of the ionising 
stars is inversely correlated with metallicity. This corre- 
lation is likely to arise from several factors, such as: (a) 
The influence of metallicity on stellar atmospheres and the 
resulting stellar spectral energy distribution, for a given 
effective temperature of the star, (b) Variations of stel- 
lar evolutionary tracks with metallicity. (c) Variations of 
IMF with metallicity. The observed spread of the correla- 
tion could result from variations in the electron temper- 
ature, starburst age and the IMF among sources in our 
sample. We discuss the causes of both the correlation and 
the scatter in this section. 

Theoretical models of hot star atmospheres indicate 
that metallicity variations affect stellar wind strengths 
and, therefore mass-loss rates, as a result of changes in 
the opacity in the line-driven winds. Changes in metallic- 
ity also induce changes in the t h e effe cts of line blocking; 
and bla nketing. IThornlev et alJ l|200fl) and IGiveon et alJ 
l)2002a|) present a discu ssion of the effects on stellar pop- 
ulations. [^^reo^^^^ find that the increase in excitation 
of Galactic H II regions with galactocentric radius cannot 
be due to this effect alone [(a) in the list above]. Increasing 
effective temperatures must strongly contribute, in addi- 



tion to the hardening of stellar SEDs for fixed T e f f . The 
same can be concluded for our gala xies by compa ring the 
observed spread in excitations to the lCiveon et alJ models. 

The effective temperature of a star-forming cluster is 
determined by at least three related factors: age, IMF, and 
metallicity. While we cannot differentiate between the in- 
dividual effects of these three parameters on our data for 
single objects, we can discuss the influence they may have 
on our abundance-excitation correlations. Stellar evolu- 
tionary tracks across the H-R diagr am also clearly displa y 
a dependence on metallicity (e.e.. ISchaller et alJ Il9921 . 
For example, lower metallicity tracks have a 'hotter' main 
sequence. This explains why high excitation regions are 
found primarily at low metallicity (as Z varies from Zq 
to Q.2Zq the ob served Ne excitation ratio increases by 
factors - 4 - 10. IThornlev et allEoOoh . Evidence for the 
Galactic Centre may suggest this effect (b) is even stronger 
(see discussion in IThornlev et~alll2000(l . 

The question of whether (c) the IMF may vary among 
clusters of differ ent metallicities remain s open. Following 
arguments as in IThornlev et alJ l)2000l) we consider such 
IMF variations a minor contributor to our correlations. 
They showed that effects of moderate IMF slope variations 
on the excitation are not large, and discussed evidence in 
favour of the most massive stars being present in the initial 
mass function of objects for a wide range of metallicities. 

We have already noted that some of our sources (par- 
ticularly the BCDs) have measured electron temperatures 
that differ from the representative value of — 5000-fC. 
However, it seems unlikely that variations in T e could 
cause a correlation between 'hardness' (effectively T e tf) 
and metallicity, as the dependen ce of abundance on e lec- 
tron temperature is weak (T e _0 7 IGiveon et alJl2002aj) . If 
T e increases with metallicity, as is suggested by the elec- 
tron temperatures we found for NGC 253, M82, II Zw 40 
and NGC 5253, then the variations in T e for our sample 
of galaxies will affect the slope of the correlation we find 
between 'hardness' and metallicity, but will not destroy it. 
In the unlikely case that T e does not systematically vary 
with metallicity, then the variations in T e in our sample 
will serve only to increase the spread in the correlation. 

The observed scatter in the correlation may also result 
from differences in both starburst ages (as argued by e.g. 
IThornlev et"aH l2000l: IGiveon et aDl2002a(l and the IMF 
for different star-forming clusters. 

5.3. Comparison to H II regions 

We also note from our 'hardness'-metallicity plots, that 
for a given metallicity, the starbursts have relatively lower 
excitation than the HII regions. Similarly low excita- 
tions have been noted i n previous studies of extragalac- 
tic star- forming regions foovon. Joseph, k. Wrightl 11994 
iDohertv et allll995t IThornlev et"al¥ 2000) . Our improved 
empiri cal database of starb ursts (this paper) and H II re- 
gions l|Giveon et alJl200 2a). reinforces this result. 
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The low excitation of galaxies has been attributed 
to low upper-mass-cu toffs (< 30Mq) and ageing effects 
l|Thornlev et alJl200fj(l . The starbursts in our sample rep- 
resent complex starburst systems with possibly multiple 
stellar populations lying within the SWS aperture, rather 
than the single stellar populations that comprise Hll re- 
gions. This difference may account for the low excitation 
since the SWS aperture may include emission not only 
from the 'youngest burst' population but also from older 
stars within the host galaxy. The net effect of this 'di- 
lution' is to lower the mean effective temperature of the 
starburst, which will result in th e starbursts having lower 
'hardn ess' ratios (cf Fig. 5 in iMartm-Hernandez et alJ 
2002bh . This lowering of the mean effective temperature 
may mirror t he effect of an ageing stellar population (as 
predicted by iThornlev et ail BoOQ 1 ) on the observed line 
ratios. 

5.4. Weak sulphur lines and the origin of low sulphur 
abundance 

We have shown in Sect. 14. 431 that the gas phase sulphur 
abundances do not follow the consistent behaviour seen in 
neon and argon. With the exception of the low metallicity 
systems, the measured S abundance is consistently smaller 
than that of Ne or Ar, by up to an order of magnitude in 
the metal-rich and dusty objects. In addition to the ques- 
tion of the origin of this deficit, the weakness of the S lines 
has practical implications for future observations of dusty 
starbursts both locally and at high redshift: sulphur lines 
are less favoured in comparison to the [Ne 11] line as tracers 
of starburst emission than one would assume from nebu- 
lar models. Instead of photoionisation models, empirical 
templates may be t he tool of choice for flux predictions. 
iGenzel et alJ lll998f) showed that low redshift ultralumi- 
nous infrared galaxies follow the trend of weak S lines 
relative to Ne, with Arp 220 being the most pronounced 
example. 

Previous studies of Galactic HII regions also 
found low S abu n dances . From infrared s pectro scopy, 
Simpson fc Rubin! Jl99dT> . ISimpson et alJ l)l995|) and 
Simpson et alJ l)l998|) derived S and Ne abundances cor- 
responding, on average, to a subsolar ratio of S/Ne, and 
suggested depletion of sulphur onto grains as a possible 
cause. This S deficienc y was recently put on a much firm er 
observational basis bv IMartm-Hernandez et all l)2002a(l 7 . 
who used ISO observations of a larger sample to dis- 
cuss Galactic abundance trends. Compared to the Ne 
and Ar abundances (with mean values similar to so- 
lar), their derived sulphur abundances are more scat- 
tered and are lower than solar by a factor ^3. For indi- 
vidual , higher-than-sola r-metallicity extragalactic H II re- 
gions, iDiaz et al. | lll99lh suggested a low S/O abundance 

7 see also Fig. El where the un der-abundance of S found by 
Martm-Herna ndez et alJ (2002a) is seen if one applies the offset 
marked by an arrow due to the differences in adopted electron 
temperature 



ratio based on optical (< 1/xm) spectroscopy, but they em- 
phasised the difficulty of determining abundances for high 
metallicity regions using optical data. These difficulties 
also apply to a possible Ar/O deficit at high metallicity 
llGarnettll20n2l . and are further reflected in the fact that 
other optical studies do not suggest a low S/O or Ar/O 
ratio for super-solar oxygen abundances ijvan Zee et alJ 
Il998l) . 

For t he low S/H abundances measured in Galactic H II 
regions, IMartm-Hernandez et al.l l|2f)02aj) propose an ex- 
planation involving uncertainties in electron temperature 
and density. While such uncertainties clearly exist, two 
arguments suggest that they are of secondary importance 
for the galaxies in our sample. First, the density measured 
from the [S III] lines is low for all galaxies in our sample 
and does not show the occasional high values that call 
the applicability of the low density approximation for de- 
termination of the S abundance in compact HII regions 
into question. Second, in addition to the inherently small 
temperature sensitivity of abundance determinations from 
fine structure lines, the remaining abundance variations 
caused by variations in temperature are expected to cor- 
relate among the three elements S, Ar, and Ne. Varying 
the adopted electron temperature from 5000 to 12500 K 
decreases the inferred abundances by factors of 2 to 2.5 
for all three elements simultaneously. Large temperature 
uncertainties therefore should not selectively lower the sul- 
phur abundance and increase the scatter in excitation ver- 
sus abundance for sulphur alone, but rather affect all three 
elements similarly. In contrast, observations single out sul- 
phur as bein g peculiar both in our data and i n the H II re- 
gion data of iMartin- Hernandez et al] l|2002ah (their Figs. 
18 and 19). 

Another issue relevant to the question of S abundance 
is the ionisation correction factor used to compute the 
total sulphur abundance. Since the S underabundance oc- 
curs in low excitation objects, the missing fraction are 
unlikely to exist in higher ionisation stages (i.e. beyond 
the 47eV ionisation potential of S 3+ ) but more realis- 
tically in S + (23eV). H o wever, the ICF calculations of 
IMartm-Hernandez et alJ l|2002al) show this effect to be 
small for the excitation range of our objects. 

We consider a direct nucleosynthesis explanation of 
the sulphur deficiency to be unlikely. Neon, argon and 
sulphur are all primary elements expected to trace each 
other reasonably well. Variations will be related to the 
different stellar mass ranges contributing to the yields for 
each element. Extreme assumptions for the initial mass 
function, and/or variations of stellar lifetimes with mass, 
can influence the ratio of neon and sulphur abundances, 
since these two elements are produced mostly in super- 
novae with high ( neon) and lower mass (sul phur) progeni- 
tors, respectively l)Wooslev fc Weaverlll995j) . However, ar- 
gon is produced over a very similar mass range as sulphur. 
Therefore, the observed abundance pattern of both neon 
and argon being 'normal', but sulphur being anomalously 
low by a large factor cannot easily be explained by differ- 
ences in the mass ranges of the SNe. 
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The calculated elemental abundances are also heav- 
ily dependent upon the atomic data and therefore er- 
rors in the published values might be responsible for 
erroneous sulphur abundances. However, we have used 
the most recently determined colli sion strengths for both 
dou bly and triply ion i sed S from ITaval fc Guotal <1999l) 
and ISaranh fc Store! l|l999ft . respectively. These papers 
use improved and more complex methods to deter- 
mine accurate collision st rengths than previous, com- 
monly used studies (e.g. iGalavis Meridoza. fc ZeinnerJ 
ll995llJohnson. Kingston, fc Duftonlll986|) . While we iden- 
tify atomic data as a possible contributor to the sulphur 
underabundance problem, the abundances we present are 
correct with respect to the most up-to-date atomic data 
available. 

A consistent explanation of both extragalactic and 
Galactic observations would be a tendency of sul- 
phur to be more strongly depleted onto dust than 
the noble gases neon and argon. However, ultravio- 
let absorption measurements in the diffuse interstel- 
lar medium often suggest most sulp hur exists in the 
gas phase llSavage fc Semb achl Il996f) with low deple- 
tion (~ 30% IWhittetl Il984t IWilms. Allen, fc McCravl 
l200f)j) . Other ultraviolet ISM studies suggest a corre- 
lation b etween sulphur dep l etion and mean hydrogen 
density (iGimdhajgkarl lHarris fc Mas Hessel Il98fit 

Ivan Steenberg fc Shulll Il988|) . with depletions reaching 
a factor of 10 for lines of sight with mean hydrogen 
densities ~ lcm~ 3 . For molecular regions, strong sul- 
phur depletion is predicted by chemical models (e. 



Dnlev. Millar, fc Williamslll980tlPrasad fc HuntresJlf 
Ru ffle et alJll999j) . Unfortunately little is known about 
the HII regions probed by our observations beyond the 
infrared studies cited above. While depletion is a natural 
explanation of the failure of the S abundance data to trace 
the noble gas abundances, a quantitative understanding 
under the conditions of our higher than solar metallic- 
ity H II regions remains to be obtained. Clearly, depletion 
would have to depend on ISM conditions to explain both 
the low metallicity dwarfs (with no S deficit both from 
our results and numerous optical studies) and the higher 
metallicity starbursts. While noting these open questions, 
we consider depletion of sulphur onto grains the best ex- 
planation for our observations. 



5.5. Wolf-Rayet Galaxies 

Wolf-Rayet galaxies exhibit signatures of Wolf-Rayet stars 
in their integrated visible spectra and have been detected 
in a variety of galaxy types (e.g. BCDs, HII galaxies, 
AGN, LINERs etc.) but not in all objec ts of a particular 
class l|Schaerer. Contini. fc Pindao1ll999l) . In general, the 
detection of Wolf-Rayet signatures in a galaxy has been 
interpreted as an indicator of young burst age (3 55 t s b Ss 
8Myr) and the presence of massive stars (M upp <; 20Mq), 
since the progenitors of Wolf-Rayet stars are postulated to 
be massive O stars. The Wolf-Rayet stage is thought to 



appear for a short time during the evolution of a simple 
stellar population. 

In our sample, Wolf-Rayet galaxies are clearly sepa- 
rated from non- Wolf-Rayet galaxies. Wolf-Rayet galax- 
ies generally have higher excitation and lower abun- 
dances than non- Wolf-Rayet objects. The current in- 
ventory of Wolf-Rayet features in external galaxies 
llSchaerer. Contini. & Pindaol ll999) is definitely not un- 
biased, but in general Wolf-Rayet signatures have been 
detected mostly in low metallicity systems. The low detec- 
tion rate of Wolf-Rayet stars in high metallicity systems is 
contrary to expectations from stellar evolutio n theory (e.g. 
lMevnetlll995|) . lLeitherer fc Heckmar] l)l995|) for example, 
model the ratio of Wolf-Rayet stars to O stars in evolving 
starbursts using the Geneva tracks, and derive an order of 
magnitude increase in the metallicity range 0.25 to 2 Zq; 
this range approximately corresponds to the metallicity 
spread in our sample. 

Recent surveys of a large nu mber of 

H II regions llBreso lin fc Kennicutl ] [2002: 

ICastellanos. Diaz, fc Terlevichl l2002t IPindao et alJ Ho02) 
located in the disks of local high metallicity galaxies have 
revealed the presence of weak Wolf-Rayet-like features 
in ~ 30% of the surveyed regions, which implies that 
Wolf-Rayet stars do in fact exist in high metallicity envi- 
ronments . 'Tliisj ; £SultM^ojTfim 

samples Jg useva. Izotov. fc Thuanl l2000t ISchaerer et all 
I2000I) and high metallicity systems incl uded in other 
Wolf-Rayet samples (e.g. NGC 3049; Vacca fc Conti 

Yet, in highly obscured star forming regions, the pres- 
ence of Wolf-Rayet stars remains to be confirmed. Our 
higher metallicity, non- Wolf-Rayet starbursts arc all well 
studied and strong Wolf-Rayet features would not have es- 
caped detection. NGC 5236 is the only exception. Despite 
its location in the excitation and abundance plots that 
is clearly consistent with th e non- Wolf-Rayet starbursts, 
iBresolin fc Kennicuttl ( 2002() have recently identified pos- 
sible Wolf-Rayet features in the optical spectrum of a nu- 
clear hot-spot (hot-spot M83-A in their paper) which lies 
within our ISO-SWS aperture. The red and blue spectra 
they present have relatively low signal-to-noise and show 
no CIVA5808 emission but do show a weak 'blue-bump' 
around 4650A which is indicative of the presence of WN 
stars. Nevertheless, it is unlikely that the blue bump is 
solely due to WN stars. NGC 5236 has a NIII A4640 / 
Hell A4686 line ratio th at is greater than unity which 
ISchmutz fc Vaccal l(l999|) have shown cannot be repro- 
duced by mixture s of known Wolf-Rayet stellar spectra. 
ISchmutz fc Vaccal suggest that a large contribution from 
Of stars could also produce such emission features. Similar 
line ratios h ave also been found i n other high metallicity 
systems (see ISchaerer et aT1l200f]j) . 

The separation of NGC 5236 with respect to II Zw 40 
and NGC 5253 is not unexpected since they represent two 
very different classes of object. NGC 5236 is large barred 
spiral galaxy which is actively forming stars on kiloparsec 
scales, whereas for II Zw 40 and NGC 5253 the majority 
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of star-formation and IR emissi on is thought to originate 
in 1 to 4 super-star clus ters fc.f. iGoriian. Turner, fc Beckl 
l200lt IBeck et alJl2002|) . Thus its location amongst the 
non-BCD galaxies is reasonable. 

There could be several reasons for the non-detection 
of Wolf-Rayet features in most of our high metallicity ob- 
jects: a special star formation history (no significant part 
of the population currently being in the brief Wolf-Rayet 
phase); dust obscuration of the regions actually hosting 
WR stars; and/or changes in spectral signatures of Wolf- 
Rayet stars at higher metallicity. 

We do not consider fine-tuned star formation histories 
to be the dominant reason for the dichotomy in our sam- 
ple. It seems very unlikely that only metal-rich objects 
are observed at extremely young or late ages (and with O 
stars still present to power the starburst). Circumventing 
this problem by postulating that the metal-rich objects are 
fully self-enriched in the current (almost terminated) burst 
would similarly require stringent assumptions on timing 
and the enrichment process. Most importantly, attempts 
using all available constraints for a detailed reconstruction 
of the star formation history in M 82, the prototype non- 
Wolf-Rayet object in our sample, predict a major part 
of the burst population to be in the 'Wolf-Ray et phase' 
<Rieke et al.lll993l iForster-Schreiher et ablEini. 

Obscuration of the regions hosting Wolf-Rayet stars is 
likely to be a significant effect. The presence of Wolf-Rayet 
stars is indicated by weak stellar emission features around 
4650A. In several of our objects, the most active star form- 
ing regions are considerably dust-obscured, and in the ab- 
sen ce of convincing near or mid-infrared Wolf-Rayet trac- 
ers ( Lumsden. Puxlev fc Dohertvlll994t) . detection of any 
'classical' optical Wolf-Rayet features is difficult. This is 
especially true if other, less obscured regions dilute the 
blue part of the visible spectrum. This is the case for M 82 
which has a stro ng post-starburst component in its opti- 
cal spectrum (e.g. iKennicuttl 1 992^ . It is worth noting that 
the high metallicity H II regions wit hin which optical Wolf- 



Rayet features have been detected ( Brcsolin & Kennicutt 
20021 ICastellanos. Diaz, fc TerlevicrJ 120021 IPindao et al " 



2002) are located within the disks of their host galax- 



ies, and may suffer less extinction (and dilution) than our 
deeply obscured infrared star forming regions. 

The role of Wolf-Rayet stars in metal rich starbursts 
also depends on still uncertain elements of the post main 
sequence evolution of massive stars. This is true for both 
the details of the evolutionary tracks, and the assignment 
of spectra to those tracks. For example, the ionising spec- 
tra of Wolf-Rayet stars may not be as hard as commonly 
assumed; this has implica tions for the spectra of compos- 
ite starburst populations llBresolin. Kennicutt. fc O arnettl 
Il999t ICrowther et alJll99flt iBresolin fc Kennicuttll2002l) . 
Additionally, evidence suggests that Wolf-Rayet stars 
in higher metallicity environments display different 
spectral signatures and weaker lines in the visi- 
ble, co mpared to those in low er metallicity environ- 
ments llMassev fc Johnson! Il998j) . Recent modelling by 
IBresolin fc Kennicuttl l)2002l) shows that Wolf-Rayet pho- 



tons have a negligible contribution to the nebular ionisa- 
tion as they are absorbed in the stellar atmospheres. The 
presence of Wolf-Rayet stars in our galaxies could be re- 
flective only of young age. Moreover, in the Galactic centre 
IThornlev et al.1 l|2000l) found that the majority of ionising 
flux originates from fairly cool (20000 to 30000 K) super- 
giants rather than the main sequence and hot Wolf-Rayet 
stars expected from a direct interpretation of the Geneva 
tracks. All this warns that a direct transfer from theo- 
retically predicted Wolf-Rayet star numbers to observed 
signatures in high metallicity environments is currently 
difficult. 

We note that the dichotomy between Wolf-Rayet 
and non- Wolf-Rayet objects in our sample could also 
be enhanced by an overestimation of the role of Wolf- 
Rayet stars in the low metallicity objects. Recent stud- 
ies suggest the existence in very young star forming re- 
gions of non- Wolf-Rayet (core hydrogen burning) stars 
that, nevertheless, show Wolf-Rayet-like broad emission 
i Massev fc Hunterll 998'l. This might additionally enhance 
the incidence of 'Wolf-Rayet features ' among the very 
youn gest, highest excitation objects l|Schmutz fc Vaccal 
Il999h . 

6. Conclusions 

We have presented mid-infrared spectral line data from a 
spectroscopic survey of a sample of starburst galaxies as 
seen by ISO-SWS. These data can be used as a reference 
database for comparison to future observations of star- 
forming galaxies with infrared telescopes such as SIRTF, 
SOFIA and Herschel. 

From this database we have investigated the excita- 
tion and abundances of the sample using fine structure 
lines of the primary nucleosynthetic products neon, argon 
and sulphur and the H-recombination lines Brackett a and 
Brackett (3. Abundance and excitation are inversely cor- 
related for Ne and Ar. In addition, a comparison to local 
H II regions shows that for a given metallicity starbursts 
are of relatively lower excitation than the H II regions. The 
excitation of mid-infrared starburst spectra is hence gov- 
erned by a combination of metallicity and other effects like 
ageing of the population, which should be accounted for 
in modelling starbursts as composite H II regions. 

An analysis of the excitation and abundance as traced 
by fine structure lines of of sulphur indicates that sulphur 
is approximately 3 times underabundant for the low exci- 
tation metal-rich galaxies relative to Ne and Ar with which 
it should be correlated since it is also a primary product 
of nucleosynthesis. We favour depletion onto dust grains 
as the most likely cause of this relative underabundance. 
The weakness of the sulphur lines may favour neon as an 
indicator of star formation in future infrared spectroscopy 
of faint galaxies. In addition the derived low sulphur abun- 
dances imply that S is not a good tracer of metallicity. 

Our sample displays a dichotomy between galaxies 
showing Wolf-Rayet features in the optical, which are 
of high excitation and low metallicity, and those with- 
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out Wolf-Rayet features, which are of low excitation and 
high metallicity. This is opposed to the expectation of 
higher Wolf-Rayet fractions at higher metallicity from stel- 
lar evolutionary models. The most plausible reasons for 
this behaviour include obscuration coupled with a lack of 
convincing mid-infrared Wolf-Rayet tracers, dilution by 
less obscured regions in the optical spectrum and changes 
in the spectral signatures of Wolf-Rayet stars at higher 
metallicity. 

The excitation probed by fine structure lines reported 
in this paper have been combined with those from a sam- 
ple of active galaxies and are presen ted in a related pa- 
per IjSturm et al . 20021 . IStiirm et alJ have constructed in- 
frared an alogues to the classical optic al excitation diag- 
nostics of lVeilleux fe Osterbrockl l)l987[) . The results show 
a clear separation between the two populations and the 
new excitation diagnostic can be used to identify pow- 
ering mechanisms of obscured sources detected in future 
infrared surveys. 
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